Population trapping due to cavity losses 
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In population trapping the occupation of a decaying quantum level keeps a constant non-zero 
value. We show that an atom-cavity system interacting with an environment characterized by a 
non-flat spectrum, in the non-Markovian limit, exhibits such a behavior, effectively realizing the 
preservation of nonclassical states against dissipation. Our results allow to understand the role of 
cavity losses in hybrid solid state systems and pave the way to the proper description of leakage in 
the recently developed cavity quantum electrodynamic systems. 
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Population trapping arises when a multilevel system 
interacting with external driving fields is frozen in a given 
state for very long times [l|, 2, 3]. The phenomenon was 
originally discovered in the dynamics of systems with at 
least three discrete levels [l[ . The same phenomenon was 
shown to appear in systems with a continuum of lev- 
els and in quantum systems interacting with quan- 
tized fields 0]. In all of these cases the mechanism is 
almost the same: population trapping occurs when there 
exists a superposition of states which decouples from the 
other states, so that its population is constant in time. In 
this case, the various transition channels corresponding 
to the states in the superposition interfere destructively, 
canceling the decay [2|. An analogous phenomenon has 
been singled out in the dissipative dynamics of an atom 
interacting with a structured reservoir. Indeed, for a two- 
level atom interacting with the quantized e.m. modes of 
a photonic-bandgap (PBG) material 0, the atomic 
population can be partially trapped in the excited state, 
when the atomic Bohr frequency is near the edge of the 
gap [1, 0, H, Q ■ In this case the trapping is due to the for- 
mation of two atom-photon dressed states, one of which, 
due to strong vacuum Rabi splitting, is protected against 
decay because the energy of the relevant transition to the 
ground state lies inside the gap, as explained in Ref. @. 
This has led to an extensive theoretical and experimen- 
tal analysis of the physics of cavities inside PBG mate- 
rials interacting either with real atoms [lfj or quantum 

dots mm. 

As discussed in Ref. [9|, what was missing until some 
years ago was the inclusion of cavity losses in the study 
of the dynamics of these systems. The cavity losses can 
be thought as due to an imperfect realization of the ex- 
perimental setup for the study of these systems. More 
precisely they arise from the imperfect alignment of the 
auxiliary waveguides which couple to the cavities accord- 
ing to the designed configurations Recently the ef- 
fect of cavity losses has been studied using a phenomeno- 
logical dissipator [l(| [Hj], identical to the one used to 
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describe the losses in cavity quantum electrodynamics 
(CQED) [l3|. However, one may wonder if and in which 
limits this phenomenological model is consistent with the 
description of a set of waveguides in the PBG material 
as a bosonic reservoir, as done in Ref. [14J. 

The scope of this paper is to provide a microscopic 
derivation of the master equation when we deal with 
CQED involving PBG materials with cavity losses, in 
the framework of an appropriate non-Markovian theory. 
We will show that the dynamics of the system shows pop- 
ulation trapping in the atomic excited state. This effect 
arises because the decay rates appearing in the micro- 
scopic dissipator are different. 

The system we study consists of a two-level atom in- 
teracting with a cavity mode, where the cavity is cou- 
pled to a bosonic environment. The interaction be- 
tween the atom and the system is described, at reso- 
nance and in units of H, by the JC Hamiltonian Hjc — 
(u)o/2)a z + ujQa^a + ft (aa+ + a^<7_) , where (a) is 
the creation (annihilation) operator of the mode, o~_ = 
\g) (e|, cr + = |e) (g\, and a z = |e)(e| - \g){g\, while \g) 
and |e) denote the atomic ground and excited state re- 
spectively [l5j]. This model is valid as long as <C ujq, so 
that one can neglect the counter-rotating terms acr_ and 
flV + . The cavity mode interacts with a bosonic reser- 
voir, with Hamiltonian Hr = uj^blpk, through the 
interaction Hamiltonian H mt — (a + a^) ^2 k gk{bk + b\), 
which has the advantage of being treatable straightfor- 
wardly in the time-convolutionlcss formalism we are ex- 
ploiting in this paper. Since the reservoir causing cav- 
ity losses is immersed in the PBG material, we expect 
its spectrum to be non-flat: to be rigorous, the master 
equation must be derived in the framework of a non- 
Markovian theory. Using the second order of the time- 
convolutionless (TCL) expansion [l6|, [TtJ , and neglecting 
the atomic spontaneous emission and the Lamb shifts, 
the master equation, in the strong atom-cavity coupling 
regime, is equal to the Markovian one under the rotating 
wave approximation (l8l . [l9j , with the important differ- 
ence that now the decay rates are time-dependent, in- 
dicating non-Markovian behavior. In this paper we fo- 
cus on the case of one initial excitation and we consider 
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a reservoir at zero temperature. In this case, the non- 
Markovian master equation for the atom-cavity system 
density operator p is: 

p(t) = -i[H JC ,p] 
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where |£? li± ) = (1/V2)(|1, .g) ± |0, e» are the eigenstates 
of Hjc with one total excitation, with energy ujq/2 ± SI, 
and |i?o) = |0, g) is the ground state, with energy —uj /2. 
The time-dependent decay rates for |J5i,_) and |J5i,+) are 
7(^0 — ^) i) and 7(^0 + fl, t) respectively. 

If the system starts from the state |0,e), i.e., if the 
atom is initially excited and the cavity is initially empty, 
from Eq. ([1]) one can obtain the system density operator 
at all times: 
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where I±(t) = J* 7(^0 ± VL,t')dt' , From Eq. © it is 
possible to compute all the populations that we will show 
in the following. Below we will study the behavior of the 
non-Markovian time-dependent rates 7(^0 ± fi, £), which, 
through the quantities I±(t), lead to population trapping. 

As a model of environment at zero temperature with 
non-flat spectrum, we consider the Lorentzian distribu- 
tion [III: 
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where a is the system-environment coupling strength, 
and A is the width of the distribution, describing also 
the inverse of the reservoir memory time. The case of 
Lorentzian spectrum is analytically treatable, while cap- 
turing important features of the non-Markovian dynam- 
ics we are interested in, i.e. the time-dependence of the 
decay rates and their different stationary values. We 
consider the case in which the spectrum is peaked on 
the frequency of the state \Ei t ~), i.e., u>i = ojq — O, 
where uj is the atomic Bohr frequency and Jl is the Rabi 
splitting due to the JC interaction. The rate 7(0;, t) for 
a generic transition with Bohr frequency uj is equal to 
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FIG. 1: Time- dependent decay rates (in units of 2fl) as a 
function of t (in units of (2S7) -1 ) for the examples we are 
considering: (i) A = 2S1/3, rates for |JSi,-) (solid line) and 
for \Ei >+ ) (dashed line); (ii) A = 2Q/V99, rates for |.Ei,_) 
(dashed-dotted line) and for (dotted line). 



7(0;, t) = 2Re {r(w, t)}, where T(uj 7 t) is related to the 
spectral density J(oj) through the relation: 
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By performing first the integral with respect to r and 
then calculating the remaining integral by means of the 
method of the residues, we obtain the following expres- 
sion for the decay rate ^(uj^t): 
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In particular, for uj\ = uj — ft and substituting ui = 
uj ± we obtain the decay rates for the two dressed 
states \Ei ±): 



7(w - tt,t) = a (l 
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for l-Ea,-) and 

7(^0 + t) 
— cos2flt 
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for \E lt+ ). 

From Eqs. ©-CO) we clearly see the general behav- 
ior of the time-dependent rates: all the rates ^{uj^t) are 
zero at t = 0, then they increase in time, till they reach 
stationary values for t S> X . These stationary values 
are proportional to J(uj), i.e., they are equal to the rates 
one obtains from a Markovian theory. For these reasons, 
the quantity A -1 can be seen as the memory time of the 
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FIG. 2: Population of the ground state of the atom-cavity 
system, as a function of t (in units of (2f2) _1 ) for the cases 
A = 2Q/3 (solid line) and A = (dashed line). The 

inset shows the short-time dynamics. 
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FIG. 3: Population of the atomic ground state, as a function 
of t (in units of (2fi) _ ), for the cases A = 2S7/3 (solid line) 
and A = 2f2/v99 (dashed line). The inset shows the very long 
time behavior. 



system-reservoir interaction and non-Markovian effects 
are expected to occur for times shorter than A -1 . 

Figure [T] shows the decay rates for two cases, corre- 
sponding respectively to A/(20) = 1/3 and A/(2ft) = 
1/V99, with a/ (2ft) = 1/10 for both cases. The con- 
dition on a assures the strong coupling regime, since 
the asymptotic decay rate of the state is 7(^0 — 

ft, 00) = a <C 2ft in both cases. As for the asymptotic de- 
cay rate of the state I Ex,+), in the first caseA/(2ft) = 1/3 
corresponds to j(ujo + ft, 00) = j(ujo — ft, oo)/10, while 
in the second case A/(2ft) = 1/V99 corresponds to 
■y(ujQ + ft, 00) = ■yiujQ — ft, oo)/100. Note that in the case 
A/(2ft) = 1/V99 the decay rate for |-Ei,+) also reaches 
negative values for short times: this is a typical feature 
of non-Markovian decay rates when the correlation time 
of the reservoir becomes large [l6| , and recently this has 
been connected to memory effects restoring coherence, 
for short times, within a quantum jump scheme [20| . 

By looking at the population P$ g of the system ground 
state |0,(7), one can monitor the loss of energy from the 
atom-cavity system. Figure [2] shows the time evolution 
of Po,g for the two cases we are analyzing. Let us first 
consider the case A/ (2ft) = 1/3 (solid line). The long- 
time dynamics is well described by a decay law which is a 
sum of two exponentials, one with rate 7(^0 — ft, 00) and 
another one with rate 7(^0 + ft, 00), slower than the first 
one. When probing the population of the ground state, 
a proper analysis of the signal should allow to point out 
this feature and to distinguish it from the purely expo- 
nential decay predicted by the phenomenological dissi- 
pator for cavity losses in strong coupling. The inset in 
Fig. [5] shows the short-time evolution of the population 
under scrutiny, where the non-Markovian effects are evi- 
dent. For very short times the rates start from zero and 
then increase linearly in time, as can be seen by expand- 
ing in power series of t the rates in Eqs. ©-0. Hence 
the short-time behavior of the population is quadratic 
in time, in agreement with the predictions of the non- 



Markovian theory [161 ] . As time increases, signatures of 
the oscillations of the rate 7(0^0 + ft, t) may appear in the 
dynamics. For A/(2ft) = 1/3 none of them is visible. For 
A/2ft = 1/V99 (dashed line in Fig. ^ the main features 
of the dynamics remain unchanged but the short-time 
dynamics shows, along with the initial quadratic behav- 
ior, also slight signatures of the oscillations in the rate, 
consisting in a non-quadratic increase in the ground state 
population, with non-monotonic derivative. Such signa- 
tures are now visible both because the memory time of 
the reservoir is much longer than in the first case and 
because the decay rate of the state |-Ei,+) also reaches 
negative values for some time intervals (see the dotted 
line in Fig. |T|) . 

In the last case, the large difference between the two 
asymptotic decay rates leads to an interesting picture of 
the long-time dynamics. Indeed the ground state popu- 
lation first reaches the value of about 50% and stays close 
to this value for a very long time, then the system starts 
to decay again and the ground state population eventu- 
ally goes to 1 (this restarting of the decay is not shown in 
Fig. [5]). It is as if the second decay channel is turned on 
after a time of the order of 7(1^0 + ft, oo) -1 . Before this 
time the decay of |£a.+) is almost completely inhibited, 
which is equivalent to the ideal situation wherein one of 
the two decay channels due to cavity losses is completely 
closed. This is at the origin of the population trapping 
due to cavity losses, as we are going to see. 

Figure [3] shows the population of the atomic ground 
state I g). This is the quantity which is usually measured 
in standard CQED experiments after the atom has left 
the cavity 13]. As in standard CQED the population 
exhibits Rabi oscillations, but the dynamics of the decay 
of the oscillations is quite different. Since the two decay 
channels have different asymptotic rates, the state 
decays before the state |J5i t +). Therefore there is a time 
interval, after the decay of |-Ei,_) and before the com- 
plete decay of \Ei t+ ), wherein the Rabi oscillations are 
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washed away by cavity losses, but some excitation is still 
present in the atom-cavity system. This is reflected in 
the fact that some population is trapped in the excited 
state of the atom for very long time, as shown in Fig. [3] 
Comparing the cases A/ (2f2) = 1/3 and A/217 = 1/V99, 
we can see that in both cases the phenomenon of trap- 
ping is clearly visible. The difference is in the amount of 
trapping and in the length of the time interval wherein 
the trapping is present. The trapping time is longer in 
the second case, since the decay of is much slower 

than in the first one. Anyway in both cases the trap- 
ping occurs for a time interval much longer than a Rabi 
period. The amount of trapping in the atomic excited 
state, in the first case, is a bit less than 20% while in 
the second case it is close to 25%. It is clear that the 
amount of population trapping increases when the ratio 
between the two rates becomes smaller. The values we 
have shown in the latter case are very close to the limit 
values for population trapping due to cavity losses. In- 
deed the limiting case is the one wherein one rate is zero. 
In this case one of the two dressed states does not decay 
and the long-time population of the system ground state 
|0, g) is 50%. The remaining population corresponds to 
the non-decaying dressed state, which is a superposition 
of the states |0, e) and |1,<?) with equal weights. This 
means that the limit value achievable for the population 
trapped in the atomic excited state is equal to 25%. 

In conclusion, we have shown that a microscopic 
derivation of the non-Markovian master equation for the 
JC model when the cavity interacts with an external en- 
vironment, allows one to predict population trapping due 
to cavity losses when the spectrum of the environment is 
not flat. It is important to note that the population trap- 
ping is a feature of the dynamics arising when the rates 



of the two dressed states are different. This means that 
even a non-Markovian extension of the phcnomenologi- 
cal model used to describe cavity losses for the JC model, 
which would contain only one single time-dependent de- 
cay rate, could not predict this effect. Only a microscopic 
derivation of the master equation for the JC model does 
allow to correctly describe cavity losses for a non-flat 
spectrum of the environment and to predict the occur- 
rence of population trapping in the dressed state, \Ei, + ). 
This state is protected against cavity losses for a very 
long time, and its photon-atom entanglement could be 
exploited, e.g., to perform experiments on the violation of 
Bell's inequality. The long lifetime could indeed allow one 
to overcome the locality loophole, making it possible to 
send the atom away from the cavity far enough to be out- 
side of the cavity lightcone [2lj . A very important point 
is that our approach paves the way to the inclusion of cav- 
ity losses in the study of more complex forms of CQED in 
PBG materials. In particular our microscopic dissipator 
could be included in some schemes of quantum feedback 
control of the quantum dynamics of the atom-cavity sys- 
tem [22j and in the study of the physics of arrays of cavi- 
ties. These latter systems have been recently exploited to 
simulate, through photon-atom dressed states, important 
models of condensed matter physics 23] . 
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